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Radiative  heat  transfer  has  a  significant  effect  on  nitric  oxide  (NO)  formation  in  turbulent  non-premixed 
flames.  Consequently,  predictive  models  of  turbulent  non-premixed  flames  must  include  an  accurate  ra¬ 
diation  submodel.  To  investigate  the  importance  of  radiation  submodels  in  modeling  NO  formation,  mul¬ 
tiscalar  measurements  of  temperature  and  species  were  coupled  with  radiation  measurements  in  a  series 
of  turbulent  non-premixed  jet  flames.  A  range  of  fuel  mixtures  were  considered  including  H2,  H2/He, 
CO/H2/N2,  CH4/H2/N2,  and  partially  premixed  CH4/air.  This  group  of  flames  represents  a  range  of 
complexity  with  regard  to  NO  formation  and  is  currently  the  subject  of  multiple  modeling  efforts.  Mea¬ 
surements  of  radiant  fraction,  temperature,  and  NO  mass  fraction  have  been  compared  with  previously 
obtained  modeling  results  for  the  H2,  H2/He,  and  CH4/air  flames.  The  results  show  that  an  emission-only 
radiation  submodel  is  adequate  for  modeling  the  hydrogen  flames  but  not  the  CH4/air  flames.  In  one 
CH4/air  flame,  the  emission-only  computations  overpredict  the  radiant  heat  loss  by  a  factor  of  2.5.  A 
comparison  of  adiabatic  and  radiative  computations  shows  that  the  inclusion  of  radiative  losses  can  reduce 
the  predicted  peak  NO  levels  by  as  much  as  57%.  An  accurate  radiation  submodel  for  hydrocarbon  flames 
must  account  for  radiative  absorption.  Spectrally  resolved  radiation  calculations  show  that  absorption  by 
COa  near  4.3  fim  is  primarily  responsible  for  the  increased  optical  density  of  the  hydrocarbon  flames.  The 
series  of  turbulent  jet  flames  considered  here  contains  a  range  of  C02  levels  and  provide  a  basis  for 
developing  a  realistic  radiation  model  that  incorporates  absorption  by  C02. 


Introduction 

The  capability  to  accurately  model  emissions  of 
nitric  oxide  (NO)  from  combustion  devices  is  an  im¬ 
portant  step  in  the  process  of  reducing  air  pollution. 
Predictive  modeling  of  NO  formation  in  turbulent 
non-premixed  flames  remains  a  significant  challenge 
because  of  the  complexity  of  the  NO  formation  pro¬ 
cess  and  the  sensitivity  of  NO  results  to  several  dif¬ 
ferent  submodels.  Attempts  to  understand  scaling 
relations  between  the  NO  emission  index  and  vari¬ 
ous  global  flame  parameters  have  had  limited  suc¬ 
cess  [1].  The  production  of  NO  depends  on  a  variety 
of  parameters  including  local  temperature,  O-atom 
concentration,  local  mixing  rates,  global  flame  resi¬ 
dence  time,  radiative  heat  loss,  and,  in  hydrocarbon 
flames,  the  fuel-rich  chemistry  of  prompt  NO  for¬ 
mation  and  reburn.  There  are  thus  multiple  sub¬ 
models  that  need  to  be  validated  before  NO  pro¬ 
duction  can  be  accurately  computed  by  turbulent 
combustion  models.  In  the  present  work,  we  con¬ 
tribute  to  the  investigation  of  the  role  of  the  radia¬ 
tion  submodel  in  computing  NO  formation. 


^Present  address:  Department  of  Mechanical  Engineer¬ 
ing,  University  of  California  at  Berkeley,  Berkeley,  CA 
94546,  USA. 


Previous  work  has  shown  that  radiation  plays  a  sig¬ 
nificant  role  in  NO  formation  in  hydrogen  jet  flames 
despite  the  relatively  low  values  of  radiant  fractions 
[2,3].  Radiation  has  also  been  shown  to  play  an  im¬ 
portant  role  in  NO  production  in  hydrocarbon 
flames  [1,4],  Radiant  fractions  in  hydrocarbon  flames 
can  be  substantially  higher  than  those  of  hydrogen 
flames  because  of  the  efficient  broadband  radiation 
from  soot.  Actual  radiant  fractions  depend  on  the 
sooting  tendency  of  the  hydrocarbon  fuel,  the  resi¬ 
dence  times  within  the  flame,  and  the  degree  to 
which  radiation  is  reabsorbed  within  the  flame. 

In  the  present  study,  we  investigate  radiative  emis¬ 
sion  and  NO  formation  in  a  series  of  non-sooting 
turbulent  jet  flames,  in  which  radiative  transfer  is 
dominated  by  gas-molecular  radiation.  The  flames 
considered  here  are  included  in  the  data  library  of 
the  International  Workshop  on  Measurement  and 
Computation  of  Turbulent  Nonpremixed  Flames 
(TNF)  [5],  which  has  been  established  to  facilitate 
collaborative  comparisons  of  measured  and  modeled 
results.  Velocity  data  and  detailed  scalar  data,  in¬ 
cluding  NO,  are  available  for  each  flame.  Fuel  com¬ 
positions  include  H2,  H2/He,  CO/H2/N2,  CH4/H2/ 
N2,  and  partially  premixed  CH4/air.  The  use  of  non¬ 
sooting  flames  eliminates  the  additional  complexity 
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TURBULENT  COMBUSTION — NOr  and  Radiation 


Fig.  1.  Measured  axial  profiles  of  the  normalized  radiant 
power  (C*)  from  four  partially  premised  turbulent  CH4/ 
air  jet  flames  (3/1  air/fuel  by  vol.)  with  different  Reynolds 
numbers.  Flame  conditions  are  given  in  Table  1.  Values  of 
C*  were  determined  using  the  total  measured  radiant 
power  for  normalization.  The  axial  coordinate  was  nor¬ 
malized  by  the  stoichiometric  flame  length. 


of  possible  soot-NO  interactions  and  reduces  inter¬ 
ference  with  the  laser-based  measurement  tech¬ 
niques.  Nevertheless,  this  series  of  fuels  corresponds 
to  an  increasing  scale  of  complexity  with  regard  to 
NO  formation.  For  example,  the  helium-diluted  hy¬ 
drogen  flames  have  low  levels  of  radiation,  such  that 
the  influence  of  the  radiation  submodel  on  NO  pre¬ 
dictions  is  minimized,  allowing  a  separate  evaluation 
of  submodels  for  coupling  turbulence  and  chemistry 
[2],  The  CO  flames  add  the  more  strongly  emitting 
and  absorbing  C02  molecule  to  the  radiation  prob¬ 
lem,  while  retaining  the  relative  simplicity  of  NO 
formation  by  the  thermal  mechanism  alone.  The 
non-sooting  methane  flames  add  complexity  in  a 
chemical  kinetic  sense,  with  the  inclusion  of  prompt 
NO  formation  and  reburn.  The  formation  of  NO  is 
strongly  affected  by  temperature  in  most  of  the  jet 
flames  considered  by  the  TNF.  Consequently,  radi¬ 
ation  submodels  must  be  validated  for  a  range  of 
flames  before  the  overall  accuracy  of  NO  predictions 
may  be  addressed. 

In  the  following  sections,  we  present  measure¬ 
ments  of  radiant  fractions  for  a  series  of  TNF  target 
flames.  For  several  of  the  flames,  computed  radiant 
fractions  were  compared  with  the  measurements.  In 
one  of  the  CH4/air  flames,  measured  spatial  profiles 
of  temperature  and  NO  mass  fraction  were  com¬ 
pared  with  both  radiative  and  adiabatic  model  cal¬ 
culations.  Radiation  calculations  using  RADCAL  [6] 
were  then  coupled  with  temperature  and  species 
measurements  to  study  the  validity  of  an  optically 
thin  assumption  in  the  radiation  submodel.  The  dis¬ 
cussion  focuses  on  the  importance  of  accurate  radi¬ 
ation  submodels  for  modeling  NO  formation  in  tur¬ 
bulent  jet  flames. 


Experimental  Methods 


Radiation  Measurements 


Radiation  measurements  were  performed  using  a 
Schmidt-Boelter-type  heat  flux  transducer  (Med- 
therm  64P-I-22)  with  a  150°  view  angle.  A  zinc  se- 
lenide  (ZnSe)  window  was  mounted  on  the  face  of 
the  radiometer  to  minimize  effects  of  convective 
cooling.  This  window  had  an  approximately  70%  op¬ 
tical  transmission  between  0.7  and  17 /an  and  passed 
the  radiation  emitted  from  the  flame,  while  isolating 
the  radiometer  from  air  currents.  The  total  radiant 
flux  emitted  by  each  flame  was  determined  following 
the  method  of  Sivathanu  and  Gore  [7],  The  radi¬ 
ometer  was  first  oriented  vertically  and  scanned 
along  a  radial  trajectory  in  the  nozzle  exit  plane  out 
to  a  distance  r  =  Lst0iCh/2,  where  Lsloich  is  the  stoi¬ 
chiometric  flame  length  determined  from  the  mea¬ 
sured  axial  profile  of  the  F avre-average  mixture  frac¬ 
tion.  The  radiometer  was  then  turned  to  face 
horizontally  toward  the  axis  of  the  jet  and  scanned 
along  the  length  of  the  flame.  The  radiant  flux  was 
then  integrated  over  a  cylindrical  surface,  which  was 
closed  at  the  base  of  the  flame  and  open  at  the  top 
to  determine  the  total  radiant  power,  Srad,  emitted 
from  the  flame. 

Examples  of  the  axial  distributions  of  radiant 
power  are  shown  in  Fig.  1  for  four  partially  premixed 
CH4/air  flames.  In  the  figure,  the  radiant  power  is 
expressed  in  terms  of  a  non-dimensional  radiant 
power,  C*,  defined  as 


C*(x/Lstoich) 


4^R3srad(x/Lstoich) 

fi'rad 


(1) 


where  R  is  the  radial  distance  from  the  burner  axis 
to  the  radiometer,  and  srad  is  the  radiated  power  as 
a  function  of  axial  position  [7],  The  C*  profiles  for 
these  four  flames  nearly  collapse  upon  one  another, 
with  the  peak  heat  release  occurring  between  %/ 
-bstoich  =  1.0  and  1.2.  These  measurements  indicate 
that  the  error  incurred  by  neglecting  radiant  loss 
through  the  top  of  the  cylinder  was  small  if  the 
height  of  the  axial  traverse  was  3  to  4  times  the  stoi¬ 
chiometric  flame  length. 

The  calibration  of  the  Medtherm  radiometer  was 
verified  using  a  copper  calorimeter.  This  verification 
ensured  that  the  factoiy  calibration,  which  was  per¬ 
formed  with  a  uniform  blackbody  source,  was  valid 
for  the  spectral  content  and  angular  distribution  of 
radiation  from  non-sooting  jet  flames.  The  calorim¬ 
eter  consisted  of  a  copper  disk  (36.1  mm  diameter 
X  3.2  mm  thickness)  with  a  known  mass  and  a  front 
surface  coating  of  flat  black  paint  having  an  absorp¬ 
tivity  of  0.93  (±0.07).  A  type-K  thermocouple  was 
attached  to  the  back  side  of  the  disk,  and  all  but  the 
front  surface  were  foam  insulated.  Radiative  and 
convective  energy  transfer  rates  were  determined 
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r/D 

Fig.  2.  Measurements  of  radial  profiles  of  ensemble  av¬ 
erage  temperature  and  species  mole  fractions  in  flame 
DCh4  at  x/D  =  45.  These  profiles  were  used  as  input  to 
the  RADCAL  calculations  for  investigating  the  optical  den¬ 
sity  of  the  flame  DCH4. 

from  measured  heating  and  cooling  curves.  The  con¬ 
vective  transfer  coefficient  was  calculated  from  the 
cooling  curves  and  used  to  correct  for  convective 
energy  loss  during  heating.  Measurements  of  radiant 
flux  were  obtained  using  the  two  probes,  and  the 
results  agreed  within  5%,  which  was  within  the  es¬ 
timated  uncertainty  of  ±  7%  for  the  calorimeter.  All 
radiant  fractions  reported  here  were  based  on  the 
original  factoiy  calibration  of  the  radiometer. 

Multiscalar  Measurements 

Experiments  to  obtain  the  multiscalar  measure¬ 
ments  included  in  the  present  paper  were  per¬ 
formed  in  the  Turbulent  Diffusion  Flame  (TDF) 
Laboratory  at  Sandia’s  Combustion  Research  Facil¬ 
ity.  Details  of  the  flow  facility  and  diagnostic  appa¬ 
ratus  have  been  described  previously  [8-11]  and  are 


not  repeated  here.  Simultaneous  point  measure¬ 
ments  of  spontaneous  Raman  scattering,  Rayleigh 
scattering,  and  laser-induced  fluorescence  (LIF) 
were  used  to  determine  temperature  and  the  con¬ 
centrations  of  CH4,  02,  N2,  H20,  COz,  H2,  CO,  OH, 
and  NO.  Fig.  2  displays  multiscalar  measurements 
of  ensemble  average  temperature  and  major  species 
for  a  radial  profile  in  a  turbulent  partially  premixed 
CH4/air  jet  flame.  Measured  profiles  such  as  these 
were  used  in  our  analysis  of  the  optical  density  of 
selected  flames.  Estimated  uncertainties  (lc)  in  av¬ 
eraged  measurements  of  scalars  used  in  the  present 
analysis  of  radiation  and  in  comparisons  with  model 
calculations  were  as  follows:  temperature  2-4%, 
C02  3-5%,  H20  3-5%,  and  NO  10-20%.  Further 
information  on  estimated  uncertainties  may  be 
found  in  the  references  given  above. 

Results  and  Discussion 

A  series  of  turbulent  axisymmetric  jet  flames  with 
a  range  of  fuel  mixtures  and  flow  conditions  were 
considered.  These  flames  are  included  in  the  TNF 
library,  and  detailed  descriptions,  including  species, 
temperature,  and  velocity  measurements,  are  avail¬ 
able  via  the  internet  [5]  and  in  the  following  refer¬ 
ences:  H2  and  H2/He  flames  [2,10,12];  CO/H2/N2 
flames  [11];  CH4/H2/N2  flames  [13];  partially  pre¬ 
mixed  piloted  CH4/air  flames  [9,14,15],  Table  1  pro¬ 
vides  a  brief  overview  of  the  flame  conditions  in¬ 
cluding  fuel  mixtures,  jet  Reynolds  numbers,  Re, 
stoichiometric  flame  length,  Lstoieh,  and  convective 
residence  times,  t.  The  convective  residence  time 
was  determined  by  r  =  Lst0ich/17jet,  where  l/jet  is  the 
jet  exit  velocity. 

Radiation  Results 

Radiant  fractions  were  measured  for  all  the  flames 
in  Table  1.  The  radiant  fraction, /rad,  is  defined  as 


TABLE  1 

Flame  Conditions  for  Turbulent  Jet  Flames 


Flame8 

Fuel  Mixture  (by  Volume) 

Re 

^-'stoich  (nun) 

t  (ms) 

ah2 

100%  h2 

10,000 

476 

1.6 

bh2 

80%  Ho.  20%  He 

9,800 

375 

1.3 

ch2 

60%  Ho.  40%  He 

8,300 

270 

1.1 

Achn 

40%  CO,  30%  H2,  30%  N2 

16,700 

197 

2.6 

Bchn 

40%  CO,  30%  H2,  30%  N2 

16,700 

340 

7.5 

adlr 

22%  CH4,  33%  H2,  45%  N2 

15,200 

512 

12.1 

Bdlr 

22%  CH4,  33%  H2,  45%  N2 

22,800 

544 

8.6 

Cch4 

25%  CH4,  75%  air 

13,400 

338 

11.5 

DCh4 

25%  CH4,  75%  air 

22,400 

338 

6.9 

ech4 

25%  CHj,  75%  air 

33,600 

338 

4.6 

fch4 

25%  CH4,  75%  air 

44,800 

338 

3.4 

fch4 

25%  CH4,  75%  air 

42,600 

338 

3.6 

'The  present  designation  of  flames  is  chosen  to  be  consistent  with  those  of  the  TNF  Workshop  [5], 
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TABLE  2 

Experimental  and  Computational  Radiant  F  ractions  for 
Turbulent  Jet  Flames 


Flame 

/rad“  ExP 

/rad  PDF 

/rad  CMC 

ah2 

9.5% 

12.3% 

11.7% 

bh2 

-5.9% 

7.3% 

6.9% 

Ch2 

-3.1% 

3.0% 

3.0% 

Achn 

3.4% 

bchn 

7.1% 

Adlr 

9.1% 

bdlr 

7.4% 

Cch4 

6.4% 

dch4 

5.1% 

10.5% 

12.5% 

Ech4 

4.1% 

fch4 

3.0% 

fch4 

3.4% 

“In  flames  BHa  and  CH.„  the  radiometer  used  to  measure 
/ra< j  had  a  sapphire  window,  which  blocked  radiation  from 
the  6.3  /rm  water  band.  The  values  have  been  scaled  based 
on  measurements  in  flame  AHa  with  the  ZnSe  window  [2], 

the  ratio  of  the  total  radiated  power,  Sra(j,  to  the 
power  released  in  the  combustion  reaction  and  is 
given  by 

-  idsh;  (2) 

where  m fuei  is  the  mass  flow  rate  of  the  fuel,  and 
/Iffcomb  is  the  heat  of  combustion.  The  measured 
radiant  fractions  are  listed  in  Table  2  for  all  the  dif¬ 
ferent  turbulent  jet  flames. 

For  the  flames  considered  here,  radiation  repre¬ 
sents  a  relatively  small  portion  of  the  total  heat  re¬ 
lease,  and  the  measured  radiant  fractions  were  less 
than  10%.  The  pure  hydrogen  flame,  flame  AHa,  had 
one  of  the  largest  radiant  fractions:  9.5%.  However, 
the  addition  of  helium  dilution  reduced  the  radiant 
fraction  to  5.9%  and  3.1%  for  20%  and  40%  dilution, 
respectively.  The  helium  dilution  decreased  the 
flame  length  and,  thus,  the  residence  time,  which 
resulted  in  a  lower  radiant  fraction.  The  reduced  sig¬ 
nificance  of  radiation  in  the  diluted  hydrogen  flames 
has  previously  been  exploited  to  examine  models  of 
turbulence-chemistry  coupling  without  the  addi¬ 
tional  complexity  of  radiation  effects  [2],  That  work 
includes  an  investigation  of  radiation  and  NO  for¬ 
mation  in  flames  Ah„-Ch2-  A  further  analysis  of  the 
radiative  properties  of  these  three  flames  is  reported 
here. 

In  the  CO/H2/N2  flames,  the  radiant  fraction  in¬ 
creased  from  3.4%  to  7.1%  when  the  nozzle  diam¬ 
eter  was  increased  by  68%  and  the  Reynolds  number 
was  kept  constant.  This  is  a  result  of  the  increase  in 
flame  volume  and  residence  time.  For  the  CH4/H2/ 


N2  flames,  the  radiant  fraction  decreased  from  9.1% 
to  7.4%  when  the  Reynolds  number  was  increased 
by  50%. 

In  the  partially  premixed  CH4/air  flames,  the  val¬ 
ues  of/ra(i  ranged  from  6.4%  for  flame  CGh4  to  3.0% 
in  flame  FCH4.  The  decrease  in  radiant  fraction  is 
primarily  a  result  of  the  reduction  in  residence  time. 
Flame  FCH4  exhibited  a  significant  amount  of  local¬ 
ized  extinction,  and  the  radiant  fraction  was  sensitive 
to  the  amount  of  extinction.  As  a  measure  of  this 
sensitivity,  we  considered  flame  F/H4,  in  which  the 
flow  velocity  was  reduced  by  5%  relative  to  that  of 
flame  FCH4.  This  slight  reduction  in  flow  velocity  de¬ 
creased  the  degree  of  extinction  and  resulted  in  a 
13%  relative  increase  in  /ra,j  from  flame  Fch4  to 
flame  F/H4.  The  sensitivity  of/ralj  to  the  initial  flow 
conditions  presents  an  additional  challenge  when 
comparing  computations  and  measurements  of 
flames  with  significant  extinction. 

Computed  Radiant  Fractions 

Calculations  and  measurements  of  radiant  fraction 
were  compared  for  the  pure  H2  and  helium-diluted 
H2  flames  and  for  flame  DCH4  of  the  CH4/air  flame 
series.  Values  of/raa  were  computed  using  both  prob¬ 
ability  density  function  (PDF)  and  conditional  mo¬ 
ment  closure  (CMC)  models  of  turbulent  jet  diffu¬ 
sion  flames  [2,16,17],  The  calculated  radiant 
fractions  are  included  in  Table  2.  The  combustion 
models  were  implemented  using  an  optically  thin  as¬ 
sumption  in  the  radiative  submodel.  This  approach 
assumes  that  the  absorption  of  radiation  within  the 
flame  is  insignificant  and  that  each  radiating  point 
has  an  unimpeded  view  of  the  cold  surroundings. 
The  flame  is  assumed  to  be  free  of  particles,  includ¬ 
ing  soot,  and  the  only  sources  of  radiant  emission  are 
gas-phase  HaO,  C02,  CO,  and  CH4. 

Radiative  transfer  in  the  H2  flames  essentially  in¬ 
volves  only  1I20.  Relatively  good  agreement  be¬ 
tween  calculated  and  measured  radiant  fractions  was 
observed  for  flames  AHo-CHo.  The  agreement  im¬ 
proved  with  increased  helium  dilution.  At  40%  he¬ 
lium  dilution,  both  the  PDF  and  CMC  calculations 
of/ralj  agreed  with  the  measurement.  This  suggests 
that  the  optically  thin  approximation  is  quite  accu¬ 
rate  for  the  diluted  hydrogen  flames  and  somewhat 
less  accurate  for  the  pure  hydrogen  flame.  A  more 
detailed  analysis  of  the  optical  density  of  flame  AHo 
is  discussed  later. 

Calculating  the  radiant  fraction  in  CH4/air  flames 
presents  an  additional  challenge  because  both  H20 
and  C02  are  major  contributors  to  the  radiative 
transfer.  In  flame  DCH4,  the  PDF’  and  CMC  calcu¬ 
lations  overpredict  /raij  by  a  factor  of  2.1  and  2.5, 
respectively.  Note  that  the  PDF  and  CMC  calcula¬ 
tions  were  performed  to  maximum  downstream  lo¬ 
cations  of  x/D  =  90  and  x/D  =  100,  respectively. 
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Axial  Position  x/D 

Fig.  3.  Comparison  of  measured  and  computed  center¬ 
line  profiles  of  temperature  in  flame  DCH4.  Ensemble  av¬ 
erage  experimental  measurements  (Exp.)  are  plotted  with 
(a)  PDF  modeling  results  using  adiabatic  (PDF  ad.)  and 
radiative  (PDF  rad.)  calculations  and  (b)  CMC  modeling 
results  using  adiabatic  (CMC  ad.)  and  radiative  (CMC  rad.) 
calculations.  An  emission-only  radiation  submodel  was 
used  in  both  the  PDF  and  CMC  radiative  calculations. 

When  the  extra  10-diameter  difference  is  consid¬ 
ered,  the  radiant  fraction  of  the  PDF  calculation  is 
actually  somewhat  closer  to  the  12.5%  value,  which 
was  obtained  for  the  CMC  calculation.  The  large 
discrepancy  between  the  calculations  and  the  mea¬ 
surements  of/ra<j  indicate  that  absorption  by  C02  is 
important  in  the  CH4/air  jet  flames. 

A  Comparative  Example 

The  incorrect  treatment  of  radiation  in  turbulent 
flame  models  can  have  a  dramatic  effect  on  the  pre¬ 
dicted  NO  levels.  As  an  example  of  this  effect,  we 
considered  modeling  results  with  and  without  radi¬ 
ation  for  flame  DCH<1.  Since  the  production  of  NO  is 
highly  dependent  on  temperature,  we  first  compared 
temperature  profiles.  In  Fig.  3,  centerline  tempera¬ 
ture  profiles  from  adiabatic  and  radiative  computa¬ 
tions  for  both  the  PDF  (Fig.  3a)  and  the  CMC  (Fig. 
3b)  models  are  displayed  along  with  ensemble  av¬ 
erage  measurements.  An  optically  thin  assumption 
was  used  in  the  radiative  calculations.  The  error  bars 
for  the  measurements  correspond  to  the  uncertain¬ 
ties  given  above.  For  x/D  <  ~40,  there  is  significant 


Fig.  4.  Comparison  of  measured  and  computed  center- 
line  profiles  of  Favre-average  NO  mass  fraction,  YNo,  in 
flame  D0H4.  Ensemble  average  experimental  measure¬ 
ments  (Exp.)  are  plotted  with  (a)  PDF  modeling  results 
using  adiabatic  (PDF  ad.)  and  radiative  (PDF  rad.)  calcu¬ 
lations  and  (b)  CMC  modeling  results  using  adiabatic 
(CMC  ad.)  and  radiative  (CMC  rad.)  calculations.  An  emis¬ 
sion-only  radiation  submodel  is  used  in  both  the  PDF  and 
CMC  radiative  calculations. 


overlap  of  the  adiabatic  and  radiative  calculations, 
and  both  PDF  and  CMC  calculations  overpredict 
the  temperature.  This  discrepancy  is  mainly  due  to 
differences  in  mixing,  since  radiation  is  not  an  im¬ 
portant  effect  for  this  portion  of  the  centerline  pro¬ 
file.  For  x/D  >  ~40,  the  radiative  calculations  pre¬ 
dict  lower  temperatures  than  do  the  adiabatic 
calculations,  indicating  the  importance  of  radiative 
heat  loss  in  this  region.  The  average  difference  be¬ 
tween  the  adiabatic  and  radiative  calculations  in  the 
upper  region  of  the  flame  is  ~190  K  and  ~130  K 
for  the  PDF  and  CMC  predictions,  respectively. 

The  difference  in  the  predicted  temperature  for 
the  radiative  and  adiabatic  calculations  has  a  signifi¬ 
cant  effect  on  the  predicted  NO  levels.  Fig.  4  shows 
a  comparison  of  adiabatic  and  radiative  calculations 
for  centerline  profiles  of  the  Favre-average  NO  mass 
fraction,  YNO-  In  Fig.  4a,  the  radiative  PDF  calcu¬ 
lation  shows  a  peak  NO  mass  fraction  that  is  37% 
lower  than  the  peak  of  the  adiabatic  calculation.  The 
agreement  of  the  radiative  PDF  calculation  and  the 
experiment  is  quite  good.  However,  the  optically 
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Fig.  5.  Results  of  a  line-of-sight  radiation  calculation  for 
a  radial  cross  section  of  the  undiluted  hydrogen  jet  flame, 
flame  AHa,  at  a  downstream  location  of  .t ID  =  34.  (a)  Spec¬ 
tral  intensity  for  emission-only  and  emission/absorption 
calculations,  (b)  Transmittance  for  the  emission/absorption 
calculation. 

thin  radiation  submodel  does  not  accurately  predict 
the  radiant  fraction  in  this  flame.  At  this  stage  of  the 
model  development,  the  agreement  between  the 
calculated  and  measured  NO  profiles  is  considered 
fortuitous  and  suggests  compensating  errors.  The 
peak  value  of  YNO  for  the  radiative  CMC  calculation 
(Fig.  4b)  is  57%  lower  than  the  peak  of  the  adiabatic 
calculation.  A  comparison  of  the  CMC  and  PDF  re¬ 
sults  in  Figs.  3  and  4  reveals  that  the  discrepancy 
between  the  adiabatic  and  radiative  predictions  of 
YNo  does  not  scale  with  the  differences  in  the  tem¬ 
perature  profiles.  This  reflects  the  complexity  of  the 
NO  formation  process  and  emphasizes  the  need  to 
independently  validate  each  submodel. 

Optical  Density  of  Flames 

The  above  results  show  the  importance  of  treating 
radiative  transport  carefully  and  suggest  that  both 
emission  and  absorption  need  to  be  considered  in 
hydrocarbon  flames.  In  order  to  better  quantify  the 


optical  density  of  these  flames,  the  RADCAL  code 
was  used  to  perform  line-of-sight  radiation  calcula¬ 
tions  for  selected  radial  measured  profiles  of  tem¬ 
perature  and  species  concentrations.  The  RADCAL 
code  computes  the  total  radiated  power  as  well  as 
the  spectrally  resolved  radiant  intensity  and  trans¬ 
mittance.  Radial  profiles  of  ensemble  average  spe¬ 
cies  concentrations  and  temperature  at  a  down¬ 
stream  location  of  x/D  =  ~Lstoic|1  were  used.  This 
location  corresponded  to  the  peak  radiative  emis¬ 
sion.  Three  representative  flames  from  the  above  se¬ 
ries  of  flames  were  considered:  flames  AHa,  DCH4, 
and  BCHN.  For  each  flame,  two  separate  RADCAL 
calculations  were  performed,  a  normal  emission/ab¬ 
sorption  calculation  and  an  emission-only  calcula¬ 
tion. 

The  radiation  calculation  for  flame  AHo  was  used 
to  investigate  the  validity  of  the  optically  thin  as¬ 
sumption  in  hydrogen  flames.  This  flame  has  the 
largest  optical  density  of  flames  AH2-CHo.  Fig.  5a 
shows  the  calculated  spectral  intensity  of  radiation 
from  H20  at  a  downstream  location  of  x/D  =  34. 
The  difference  between  the  optically  thin  and  opti¬ 
cally  thick  calculations  is  relatively  small.  The  cal¬ 
culated  value  of  total  radiated  power  for  the  emis¬ 
sion-only  calculation  is  only  13%  higher  than  that  of 
the  emission/absorption  case.  Fig.  5b  shows  that  the 
small  amount  of  absorption  by  H20  resulted  in  a 
minimum  transmittance  of  74%  at  6.7  //m.  These 
results  indicate  that  the  optically  thin  assumption  is 
adequate  for  the  hydrogen  flames. 

The  measured  species  and  temperature  profiles 
used  in  the  RADCAL  calculation  for  flame  DCH4  are 
shown  in  Fig.  2.  The  calculated  value  of  total  radi¬ 
ated  power  for  the  emission-only  case  is  39%  higher 
than  that  of  the  emission/absorption  computation. 
This  indicates  that  optical  absorption  is  important  in 
flame  DCH<1.  The  calculated  spectral  distributions  of 
the  emission  intensity  and  transmittance  for  flame 
DCHj  are  plotted  in  Fig.  6.  Results  for  the  4.3  pm 
band  of  C02  show  a  significant  effect  of  absorption 
over  the  measured  flame  profile.  Fig.  6b  shows  that 
the  transmittance  for  the  emission/absorption  cal¬ 
culation  drops  to  a  minimum  of  17%.  These  results 
indicate  that  significant  errors  in  the  calculated  ra¬ 
diant  fraction  should  be  expected  if  an  optically  thin, 
or  emission-only,  assumption  is  used  in  treating  C02 
radiation  from  these  CH4/air  flames. 

The  CO/H2  flames  generate  C02  levels  approxi¬ 
mately  twice  those  in  the  CH4/air  flame.  Conse¬ 
quently,  absorption  by  COz  is  expected  to  be  even 
more  significant.  A  radiation  calculation  was  per¬ 
formed  for  flame  BCHN,  which  is  the  larger  of  the 
two  CO/LL  flames.  The  results  of  the  RADCAL  cal¬ 
culation  are  shown  in  Fig.  7.  As  was  the  case  with 
flame  DCH4,  the  most  prominent  feature  is  the  C02 
band  near  4.3  pm.  The  peak  spectral  intensity  in  the 
emission/absorption  calculation  is  only  13%  higher 
than  that  for  flame  DCH4  (see  Fig.  6a).  For  the 
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Fig.  6.  Results  of  a  line-of-sight  radiation  calculation  for 
a  radial  cross  section  of  flame  DCH4  at  a  downstream  lo¬ 
cation  of  x/D  =  45.  The  computation  was  performed  using 
the  measured  radial  profiles  of  species  and  temperature 
shown  in  Fig.  2.  (a)  Spectral  intensity  for  emission-only  and 
emission/absorption  calculations,  (b)  Transmittance  forthe 
emission/absorption  calculation. 

CO/H2  flame,  however,  the  difference  between  the 
emission-only  and  emission/absorption  calculations 
is  much  more  significant.  The  total  emitted  intensity 
for  the  emission-only  calculation  is  2.2  times  that  of 
the  emission/absorption  case.  Fig.  7b  shows  that  the 
transmittance  associated  with  4.3  jum  band  of  C02 
has  a  minimum  of  4.5%.  Clearly,  the  computation  of 
NO  formation  in  the  CO/H2/N2  flames  must  in¬ 
clude  the  radiative  absorption  of  COa. 


Conclusions 

The  importance  of  radiation  in  modeling  N O  for¬ 
mation  was  investigated  for  a  series  of  turbulent  non- 
premixed  jet  flames  using  combined  multiscalar 
diagnostics  and  radiometric  measurements.  Mea¬ 
surements  of  radiant  fraction,  temperature,  and  NO 
mass  fraction  were  presented  and  compared  with 
previous  modeling  results.  The  models  included  an 


Fig.  7.  Results  of  a  line-of-sight  radiation  calculation  for 
a  radial  cross  section  of  flame  BChn  of  the  CO/H2/N2 
flames  at  a  downstream  location  of  x/D  =  50.  (a)  Spectral 
intensity  for  emission-only  and  emission/absorption  cal¬ 
culations.  (b)  Transmittance  for  the  emission/absorption 
calculation. 


emission-only  radiation  submodel,  which  we  found 
to  be  adequate  for  pure  hydrogen  flames  and  he¬ 
lium-diluted  hydrogen  flames.  However,  in  hydro¬ 
carbon  flames,  the  optically  thin  assumption  was  in¬ 
appropriate.  In  a  partially  premixed  CH4/air  flame, 
computations  overpredict  the  radiant  fraction  by  a 
factor  of  2.5.  A  comparison  of  results  from  radiative 
and  adiabatic  computations  in  this  same  flame  dem¬ 
onstrated  the  sensitivity  of  the  predicted  NO  levels 
to  radiation.  The  inclusion  of  radiative  losses  with  a 
radiant  fraction  of  12.5%  resulted  in  a  57%  reduc¬ 
tion  in  predicted  NO  levels. 

Further  analysis  of  the  spectral  characteristics  of 
the  radiant  emission  and  transmittance  confirmed 
that  an  optically  thin  assumption  is  inappropriate  for 
both  the  CO/H2  and  partially  premixed  CH4/air  jet 
flames.  The  primary  source  of  radiative  absorption 
was  C02  in  the  4.3  jum  band.  The  need  to  incorpo¬ 
rate  optical  absorption  into  radiation  submodels 
complicates  modeling  efforts.  A  full  radiative  trans¬ 
fer  calculation  is  computationally  expensive,  and  an 
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alternative  method  of  including  C02  absorption  is 
needed  for  hydrocarbon  jet  flames.  Previous  efforts 
to  incorporate  radiation  into  modeling  of  turbulent 
nonpremixed  flames  include  work  with  strongly  ra¬ 
diating  acetylene/air  flames  by  Gore  et  al.  [18].  The 
measurements  and  analysis  presented  here  will  serve 
as  a  basis  for  developing  a  valid  radiation  submodel 
for  the  TNF  library  of  flames.  Once  the  radiation 
submodel  has  been  appropriately  addressed,  it  will 
be  feasible  to  study  the  effects  of  various  aspects  of 
the  turbulence-chemistry  models  on  NO  formation. 
Since  NO  formation  is  less  complex  in  CO/H2 
flames  than  in  CH4/air  flames,  it  may  be  useful  to 
first  address  the  submodels  for  CO/H2  flames. 
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COMMENTS 


Stephen  B.  Pope,  Cornell  University,  USA.  The  absorp¬ 
tion  has  been  estimated  based  on  mean  profiles,  that  is, 
with  the  neglect  of  turbulent  fluctuations.  Can  you  quantify 
the  impact  of  turbulent  fluctuations  on  your  conclusions. 

Author’s  Reply.  The  use  of  ensemble  average  radial  pro¬ 
files  of  temperature  and  species  concentrations  in  the 
RADCAL  calculations  did  not  account  for  the  turbulence/ 
radiation  interactions.  The  quantitative  evaluation  of  tur¬ 
bulence/radiation  interactions  requires  measurements  of 
the  instantaneous  radial  profiles  of  temperature  and  species 
concentrations.  The  capability  to  perform  such  measure¬ 
ments  is  currently  under  development  [1]. 

Previous  results  comparing  predictions  of  mean  and  sto¬ 
chastic  properties  in  both  hydrogen/air  and  methane/air 


diffusion  flames  indicate  that  the  turbulence/radiation  in¬ 
teractions  are  important  to  consider  but  do  not  alter  the 
conclusions  of  our  analysis  of  the  optical  density  of  turbu¬ 
lent  jet  flames  [2,3] . 
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